Multidimensional optical fractionation with holographic verification 
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The trajectories of colloidal particles driven through a periodic potential energy landscape can 
become kinetically locked in to directions dictated by the landscape's symmetries. When the land- 
scape is realized with forces exerted by a structured light field, the path a given particle follows 
has been predicted to depend exquisitely sensitively on such properties as the particle's size and 
refractive index These predictions, however, have not been tested experimentally. Here, we describe 
measurements of colloidal silica spheres' transport through arrays of holographic optical traps that 
use holographic video microscopy to track individual spheres' motions in three dimensions and 
simultaneously to measure each sphere's radius and refractive index with part-per-thousand resolu- 
tion. These measurements confirm previously untested predictions for the threshold of kinetically 
locked-in transport, and demonstrate the ability of optical fractionation to sort colloidal spheres 
with part-per-thousand resolution on multiple characteristics simultaneously. 



Optical fractionation is akin to sifting particulate mat- 
ter with a sieve made of light [l|,S,[3. Colloidal particles 
flowing with or through a viscous medium such as water 
are deflected into different directions by their differing 
interactions with a structured light field [H . Unlike pho- 
tophoresis in which particles are deflected down the 
optical axis by radiation pressure, optical fractionation 
relies j3n_ optical forces directed transverse to the optical 
axis ^5,, 6] and on the symmetries and periodicities of the 
structured light field to select different fractions from an 
initially heterogeneous sample 0, H, S B H 0; U I • 

Pioneering studies 0, i i a i Hlj l have suggested 

that optical fractionation may offer unparalleled selectiv- 
ity for particle size and refractive index. Lacking suffi- 
ciently precise methods to test these predictions directly, 
experimental studies either have focused on sorting poly- 
disperse samples of particles with easily distinguishable 
properties 

0,1,111 

or else have used indirect approaches 




13|. 



to assess the method's finesse 

This Letter describes optical fractionation measure- 
ments on nominally monodisperse colloidal samples 
whose performance is verified directly using the high- 
resolution single-particle characterization and tracking 
capabilities [ij] of holographic video microscopy [l5i] . 
These measurements demonstrate optical fractionation 
on the basis of size and refractive index with part-per- 
thousand resolution and agree quantitatively with predic- 
tions derived from limiting arguments for colloidal trans- 
port through optical trap arrays 0, Hi- They thus cast 
new light on the robustness of kinetically locked-in trans- 
port against the inevitable disorder in practical imple- 
mentations. 

Our experimental system is depicted schematically in 
Fig. [TJ Aqueous dispersions of colloidal silica spheres 
are driven by a syringe pump (Harvard Apparatus 
PHD 2000) at a constant flow rate of 4 //1/h through 
a rectangular microfluidic channel formed by bonding 
a poly(dimethylsiloxane) casting to a glass microscope 
cover slip The resulting channel is 3 cm, 30 /im 



holographic 
trapping 



FIG. 1: (color online)(a) Combined holographic optical trap- 
ping and holographic video microscopy system for tracking 
and characterizing colloidal spheres during optical fraction- 
ation. The objective lens focuses computer-generated holo- 
grams projected with a spatial light modulator (SLM) into 
the sample. It also collects holographic images of the sam- 
ple that are recorded with a video camera. Trapping and 
imaging light are separated with a dichroic mirror, (b) Image 
of 103 holographic optical tweezers in the microscope's focal 
plane. The traps are arranged in two domains of oblique lat- 
tice points with a = 3.27 /im and h — 2.025 ^m. One unit 
cell is highlighted. The domains are oppositely inclined so 
that 6i[io] = ±30°, 6'[oi] = t38° and 9[ii] = t13° with re- 
spect to the flow direction. Dashed boxes indicate input and 
output regions. Traces show trajectories for four representa- 
tive particles measured at 1/30 s intervals, (c) Normalized 
holographic image of 1.5 ^m diameter silica spheres in the 
same field of view. Inset: Fit to the Lorenz-Mie theory of the 
outlined sphere's holographic image. 



deep and 1 mm wide. Spheres travel along the channel 
in a parabolic Poiseuille flow profile, their speed depend- 
ing on their height in the channel. 

The flowing spheres interact with an array of holo- 
graphic optical traps 17, [l^l whose focused intensity 
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pattern is recorded in Fig. [Ijb). The pattern of 103 
optical tweezers Q is created from a single laser beam 
by imprinting a computer-generated hologram [3| onto 
its wavefronts using a phase-only spatial light modulator 
(SLM; Hamamatsu X8267-16). The hologram is relayed 
to the objective lens (Nikon Plan-Apo, lOOx, NA 1.4, oil 
immersion) of an inverted optical microscope (Nikon TE- 
2000U), which focuses it into diffraction- limited traps at 
the midplane of the sample [l^. Figure [T](b) was created 
by placing a front-surface mirror in the microscope's fo- 
cal plane, collecting the reflected light with the objective 
lens, and recording it with the microscope's video cam- 
era (NEC TT324A II). The traps are powered by a sohd 
state laser (Coherent Verdi 5W) operating at a vacuum 
wavelength of A = 532 nm. 

The optical tweezer array was designed to sort an input 
distribution of colloidal spheres into two fractions that 
are kinetically locked in to two distinct lattice directions, 
and thus are deflected by the array in opposite directions. 
The array consists of two oblique lattice domains whose 
[10] axes are oppositely inchned at 6[io] — ±30° with re- 
spect to the flow direction. Photometry performed on 
the traps' images suggests that their intensities vary by 
roughly 10%. This is consistent with calibration mea- 
surements based on the thermally-driven fluctuations of 
spheres trapped in each tweezer under zero-flow condi- 
tions Trajectories that become kinetically locked 
in to the [10] direction, such as the two measured exam- 
ples denoted by circles in Fig.[TJb), are deflected toward 
the apex of the array before flowing downstream. Other 
particles that interact too weakly with the traps to be 
deflected by 30° may still become locked in to the [11] 
direction, which is inclined at d[ii] — Tl3° in the two 
domains. Particles traveling along this direction are de- 
flected away from the apex to the sides of the array, as in- 
dicated by the two typical trajectories denoted by squares 
in Fig.lHb). The static array of traps therefore sorts the 
sample into two distinct fractions, one concentrated at 
the apex of the array, and the other deflected to its sides. 

The objective lens also is used to holographically 
track and characterize individual spheres as they travel 
through the array. To do this, the sample volume is 
illuminated with the coUimated 3 mm diameter beam 
from a HeNe laser operating at a vacuum wavelength 
of 632.8 nm. The laser's irradiance is on the order of 
1 nW/^m^, which is less than 10~* of an individual trap's 
irradiance and is comparable to that of a conventional in- 
candescent illuminator. Light scattered by a sphere inter- 
feres with the unscattered portion of the beam to form an 
interference pattern in the microscope's focal plane. This 
is magnified and its intensity recorded by the video cam- 
era with a calibrated spatial resolution of 91 nm/pixel 
at 30 frames/s. Figure [Tfc) shows a detail of a typical 
holographic snapshot of 1.5 fim diameter silica spheres 
normalized by a background image to eliminate speckle 
and other static nommiformities |14i. 



Each snapshot is analyzed ij] using the Lorenz-Mie 
theory of light scattering [2l| to obtain each colloidal 
sphere's three-dimensional position, rp{t), relative to the 
center of the focal plane at time t, its radius Op, and 
its complex index of refraction, rip. Rigorous error esti- 



mates jl4lj yield nanometer resolution for Vp and Op, and 
part-per-thousand resolution for Up. Such precise single- 
particle characterization makes possible a direct test of 
optical fractionation's resolution. 

A particle becomes kinetically locked-in to one of the 
array's crystallographic directions because its interaction 
with one trap defiects it toward the next trap along that 
direction, and so on through the lattice P. Isl. Fzl [^. fill. [2^ . 
Modeling an individual trap as a Gaussian potential well 
of width (Tp and depth Vp separated from its neighbors by 
a distance &, the maximum angle ^max by which a particle 
driven by force Fp can be deflected has been estimated 
d 0, [H to satisfy 



< sin 
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^ ■ (1) 



This is a general result that applies equally well to ar- 
rays of potential barriers as to potential wells 0, [l^] ■ In 
the particular case of a micrometer-scale colloidal sphere 
interacting with optical tweezers, the effective well depth 
depends on both the particle's size and refractive index, 
and is given approximately by 0, [l^ , 



Vp^ 



Up + 2n^ 



(2) 



where rim is the refractive index of the medium and P is 
the laser power in an individual trap. The well's effective 
width 0, 0], ap = ap + A^/(4n^), also depends on the 
particle's size because larger spheres encounter a trap's 
localized light field from greater distances. Finally, the 
driving force due to Stokes drag in a fiuid of viscosity 77 
flowing at velocity v is 



Fp — Qirrjap v, 



(3) 



assuming that hydrodynamic coupling to boundaries and 
other spheres is weak enough to ignore 24| . 

Equations ([T]), ^ and ^ define conditions for 
marginally locked-in transport that divide the {ap^rip) 
plane into a class of particles that can be deflected to 
angle 9 and a complementary class that cannot. This 
division falls along the curve 
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FIG. 2: (color online) High-resolution optical fractionation of mono disperse colloidal silica spheres demonstrated in the measured 
probability density p{ap,np) of particles with radius ap and refractive index Hp. The input sample, whose distribution is shown 
in (a), is separated into a smaller, low-index fraction in the side regions (b) and a larger, high-index fraction at the apex (c). 
The dashed curves show the marginally locked-in condition nc(op) for the array's [10] direction. Plotted squares and circles 
indicate the properties of the four representative spheres whose trajectories are plotted in Fig. [TJb). 



depends on particle size through the ratio Op/b and 
B = — — b sm d 



(6) 



depends only on control parameters. Particles satisfying 
Up > nc{ap) should become locked in to the lattice di- 
rection at angle 9, and those with Up < nc{ap) should 
escape, perhaps to become locked in along other lattice 
directions at smaller angles. This result provides a basis 
to tune the selectivity of optical fractionation for sorting 
by size and refractive index. 

The data in Fig. [2] show results obtained with a sin- 
gle sample of monodisperse silica spheres (PolySciences 
Catalog #24327, Lot #600424). Lorenz-Mie characteri- 
zation of 6000 randomly chosen spheres yields a sample- 
averaged radius of Up = 0.715 ± 0.021 /zm, which is sig- 
nificantly smaller and more monodisperse than the man- 
ufacturer's specification of 0.75 ± 0.04 ^m. The sample- 
averaged refractive index of Up = 1.418 ± 0.004 is consis- 
tent with typical values for colloidal silica spheres [l3 |. 
The distribution of input particle properties in Fig. ^a) 
was computed from single-particle results using an opti- 
mal non-parametric density estimator (25l |. It reveals an 
anticorrelation between particle size and refractive index 
that is typical for porous colloidal spheres grown by emul- 
sion polymerization and is not observed in fully compact 
particles such as fluid droplets [l3]- 

The homogeneous input sample was flowed at v — 
34 ± 2 nm/s through the array of traps. This is slow 
enough that blurring during the camera's 1 ms expo- 
sure time does not de gra de holographic particle tracking 
and characterization [ij]. Although particles approach 
the array at a range of heights and thus at a variety of 
speeds, they are drawn to the channel's midplane by op- 
tical forces, and thus pass through at a uniform speed. 



This vertical flow focusing is made apparent by the three- 
dimensional tracking capabilities of holographic video mi- 
croscopy. The concentration of particles is low enough 
that no more than eight occupy the array at any time. 
Collisions are rare at this low occupancy and no effort 
was made to exclude their effects from our results. 

The trap array's [10] lattice spacing of 6 = 2.03 /im 
was chosen so that roughly half of the sample would be 
kinetically locked-in according to Eq. (j4|) with a laser 
power of P = 3 mW/trap. The critical condition for 
locked-in transport along the [10] lattice direction pre- 
dicted by Eq. ^ is plotted as a dashed curve in Fig. [2] 
with no adjustable parameters. The [01] lattice spac- 
ing of a = 3.27 ^m was selected to be large enough to 
prevent locked-in particles from escaping [7|, but small 
enough that the low-angle [11] direction would be able 
to entrain and deflect all of the particles in the sample. 
The array thus preferentially deflects particles toward its 
apex, but provides an alternative route with opposite de- 
flection for particles that interact too weakly to lock in 
to its principal axis. The two fractions thus are spatially 
separated from each other and also from those particles 
that bypass the array's catchment basin altogether. 

The characteristics of particles flowing into and out of 
the array were assessed within the rectangular regions in- 
dicated in Fig.lHb), with results from the two side regions 
being combined. Particle-by-particle characterization in 
the side and apex regions, plotted in Figs. [^Jb) and (c), 
reveals clear separation into two distinct and sharply de- 
fined sub-populations, the larger, higher-index particles 
being deflected to the apex and the smaller, lower-index 
particles being deflected to the sides, in excellent quan- 
titative agreement with Eq. ([?]). The mean radius of 
the focused population is slightly larger than that of the 
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other fraction (0.716±0.021 fj.m versus 0.710±0.027 /im). 
A far more substantial distinction is observed in the re- 
fractive indexes of the sorted populations. The focused 
spheres have a mean refractive index of 1.420 ± 0.003 
whereas the dispersed fraction has a mean refractive in- 
dex of 1.416 ± 0.003. The properties of the two spatially 
separated populations differ with better than 99.99 per- 
cent confidence according to the Wilcoxon rank-sum test. 

The results in Fig. [2] demonstrate that optical frac- 
tionation is inherently a multidimensional sorting tech- 
nique whose selectivity for size and refractive index 
can be dynamically optimized by tuning the proper- 
ties of the optical trapping pattern. The present re- 
alization demonstrates discrimination on the basis of 
refractive index with part-per-thousand resolution un- 
der the set of conditions that we considered. Similarly 
good results, including agreement with Eq. ([4]), were ob- 
tained for a monodisperse sample of micrometer-diameter 
polystyrene spheres. These observations provides exper- 
imental support for the suggestion 0, [ill, IS] that the 
large number of particle-trap interactions in a typical 
trajectory renders kinetic lock-in robust against the noise 
and quenched disorder that arise in practical implemen- 
tations. A rigorous explanation for this robustness re- 
mains an outstanding challenge, as does the quantitative 
predictive power of Eq. ([4]), which was derived with semi- 
quantitative limiting arguments 0, 0] . 

The inhomogeneous trapping pattern used for this 
study implements a prismatic mode of optical fraction- 
ation P, |Tl[ in which different fractions become kineti- 
cally locked in to different crystallographic directions and 
thus are sorted into spatially separated streams. More 
complex trapping patterns should enable still more so- 
phisticated multi-dimensional sorting with similarly large 
displacements for finely resolved populations. The very 
high-resolution fractionation by refractive index that we 
have demonstrated should have immediate applications 
in photonics, where colloidal particles' optical properties 
increasingly are being exploited, and also in medical di- 
agnostics, where pre-selecting probe particles for precise 
optical properties may reduce the time, effort and ex- 
pense of bead-based molecular binding assays [1J| . The 
throughput for such applications can be increased by pro- 
portionately increasing the flow speed and laser power. 

This work was supported by the MRSEC program of 
the National Science Foundation through Grant Number 
DMR-0820341. 
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